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Abstract: We explore an extension of the Standard Model by an additional 17(1) gauge 
group and a complex scalar Higgs portal. As the scalar is charged under this gauge factor 
this simple model supplies a vector dark matter candidate satisfying the observed relic 
abundance and limits from direct dark matter searches. An additional Higgs-like state, 
that may be heavier or lighter than the observed Higgs, is present and satisfies LEP and 
LHC bounds whilst allowing for absolute stability of the electroweak vacuum in a range of 
the parameter space. 
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1 Introduction 

The discovery of a 125 GeV scalar [1, 2] completes the Standard Model (SM) of particle 
physics and appears to confirm some rather intriguing features of nature. Firstly, the Higgs 
potential of the SM develops an instability at high field values, [3-10]. Secondly, the theory 
on its own does not account for nearly 95% of the matter in the universe, known as dark 
matter (DM) [11] and dark energy. 

In this paper we explore a simple extension of the SM by an extra U{l)x gauge 
symmetry factor, supplemented with an additional complex scalar charged solely under 
this U{l)x- It was studied by [12-15] and by [16-18] in its non-Abelian version (for other 
realizations of vector fields in the context of dark matter see [19, 20]). The model can (if 
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the real component of the complex scalar obtains a vev) generate a viable massive vector 
dark matter candidate. It also provides an additional Higgs mass eigenstate and introduces 
additional freedom in the theory to completely alleviate the issue of vacuum stability. 

In this paper we use a collection of HEP-tools and software to give an accurate account 
of the currently viable parameter space of such a model. In particular we combine SARAH 
[21-24], Pyrate [25] and SPheno [26, 27] to explore stability at 2-loops and micrOMEGAs 
[28-31] for a detailed study of vector dark matter for this model. We also compare our 
results to various experiments including LUX [32] and XENONIOO [33]. 

The paper is organized as follows. In section 2 the model is defined and a detailed 
exposition of the parameters of the model is presented. In section 3 the renormalisation 
group equations of the theory are used to constrain the parameter space from physical 
consideration such as positivity, stability and the absence of Landau poles in the theory. 
In section 4 collider constraints are applied and the dark matter abundance constraint is 
imposed on the parameter space that remained. A detailed study of the DM-nucleon cross 
section is also performed. Einally in section 5 the findings are summarized. 

In Appendix A we give a detailed account of the tree-level calculation of the relic 
abundance. Appendix B supplies the SARAH model files which fixes the conventions and 
allows for the generation of code for RGES and for the code used in micrOMEGAs [28-31]. 

2 A model of vector boson dark matter 

The vector dark matter (VDM) model [12-16] is an extension of the SM by an addi¬ 
tional U{l)x gauge symmetry together with a complex scalar field S', whose vev generates 
a mass for this 1/(1)’s vector field. The quantum numbers of the scalar field are 

S= ( 0 , 1 , 1 , 1 ) under t/(l)Y x SU(2 )l x St/(3)c x ( 2 . 1 ) 

None of the SM fields are charged under the extra gauge group. In order to ensure stability 
of the new vector boson a Z 2 symmetry is assumed to forbid t/(l)-kinetic mixing between 
U{l)x and U{1)y- The extra gauge boson and the scalar S field transform under Z 2 
as follows 


—>• — , S ^ S*, where S = 4)e^^, so (p ^ 4>, u —)• —a. (2.2) 


All other fields are neutral under the Z 2 . 

The charge neutral vector bosons of U{1)y & SU{2)i mix and the model 

leads to ( 7 , Z). The mixing matrix is taken to be 


B \ _ ( cos 9w — sin \ / 7 

IT3 / i singly cosily ] \ Z 


where 


1 


cos 6w = 




, singly = 


9'^ +9'^ 


^9' = ^3/551 and 92. 


(2.3) 


(2.4) 
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At leading order the vector bosons masses are given by; 


Mw = ^gv, Mz = ^Vg^ + and Mz' = gxVx, (2.5) 

where v and Vx are H and S vacuum expectation values (vev’s): {{H),{S)) = -^{v,Vx)- 
The scalar potential for this theory is given by 

y = + \h\H\^ - 4|5|2 + As|5|4 + k\S\^\H\\ (2.6) 

It will also be useful to define, for future reference, the parameter \sm = = 0.13, 

where Mh = 125.7 GeV. 

The requirement of positivity for the potential implies the following constraints that 
we impose in all further discussions: 

Xh > 0, Xs > 0, K > -2 ^JXhXs- (2.7) 


Hereafter the above conditions will be referred to as the positivity or stability conditions. 

It is easy to find the minimization conditions for scalar fields (without losing generality 
one can assume v,Vx > 0 ): 

{2Xhv^ + — 2fi\)v = 0 and {kv^ + 2Xsv^ — 2fig)vx = 0 ( 2 . 8 ) 


If fijj 5 < 0 the global minimum at (0, 0) is the only extremum. For 5 > 0 the point (0,0) 
is a local maximum of the potential, in this case ( 0 , and ( ^^)^ , 0 ) are global minima 
if > AXhXs, otherwise they are saddle points and the global minima are determined by 


2 _ 2Kn‘g 2 _ 2/t/i|^ 

4A//A5 — ^ AXjjXg — K? 


(2.9) 


For the VDM model only the latter case is relevant, since both vevs need to be non-zero 
to give rise to the masses of the Standard Model fields and dark vector boson. Both scalar 
fields can be expanded around corresponding vev’s as follows 


S 


1 

71 


{vx + (t>s + i(^s) 




1 

71 


{v + cpH + icrn) 


where H 



( 2 . 10 ) 


The mass squared matrix for the fluctuations {4>h,4’s) and their eigenvalues read 


= 


2Xhv^ kvVx 
KVVx 2Xsvl 


Ml = 


Xhv^ + ^‘s^x ~ ‘2 ^XhXsv‘^v‘ 1 + X^^v'^ + 


( 2 . 11 ) 


The matrix could be diagonalized by the orthogonal rotation R, such that = 

R~^Ai^R. The convention adopted for the ordering of the eigenvalues and for mixing angle 
a is the following 




(K “ ) 
I » Kj 


R = 


( cos a — sin a 
sin a cos a 



( 2 . 12 ) 
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where = 125.7 GeV is the mass of the observed Higgs particle. Then we obtain 


sin2a = ^ sign(AsM - Ah)(A<;, - 

V(X'fi-Xi2)2 + 4(Xa2 v'(>l'fl-Xi2)2 + 4(Xa2 

Note that since vev of H, if fixed at 246.22 GeV, with k = 0 (no mass mixing) and 
7 ^ ^SM it is only /12 which can have the observed Higgs mass of 125.7 GeV. Even 
though the mass matrix is diagonal in this case, however in order to satisfy our convention 
that = 125.7 GeV a rotation by a = ± 7 r /2 is required in such a case. 

There are 5 real parameters in the potential: ^h-, fJ-s, A 5 and k. Adopting the 
minimization conditions (2.8) fj-n, IJ-s could be replaced by v and Vx- The SM vev will be 
fixed at u = 246.22 GeV. Using the condition Mh^ = 125.7 GeV, v‘1 could be eliminated 
via ( 2 . 11 ) in terms of Xh, k, Xs, Xsm = 


2 _ 2 4AgM(AH - Xsm) 
4:Xs(Xh - Xsm) - k? 


(2.14) 


Therefore eventually there are 4 independent unknown parameters in the model {Xh, k, Xs, gx), 
where Qx is the U{l)x coupling constant. 

It is important to notice that positivity of implies for Xh > Xsm that 


Air > VT—h Asm (2-15) 

4As 

Applying the fact that = 2{Xhv‘^ + Agu^) together with (2.14) one finds the 

following universal formula for the mass of the non-standard Higgs: 


^2 _ 2 2(Ag - Xsm){4^XhXs - K^) 
4As(A_h- - Xsm) - 


(2.16) 


It is easy to see that positivity of is guaranteed if the following conditions are satisfied: 

2 

• for Xh < Xsm, Xh > 

2 

• for Xh > Xsm, > ^ + Xsm (same as (2.15)). 

From (2.11) one can easily derive the following useful inequalities 

M+ ^ Xhv^ + Xsv^ + \Xhv^ — Xsv^\ = Max(2Aisu^,2Asu^), (2-17) 

< Xhv^ + Xsvl — \Xhv‘^ — Xsv^l = Min( 2 AHf^, 2Xsvl). (2-18) 

Therefore for Xh > Xsm we have M+ > M/j^, so the lighter scalar is the SM-like Higgs 

particle, while when Xh < Xsm the observed Higgs particle is the heavier one. Con¬ 

sequently Xh is the sole parameter that distinguishes between these two scenarios. To 
illustrate that behavior we show in fig. 1 contours of non-standard Higgs masses in the 
plane {k,Xh)- Note the presence of = 0 contour (the most external blue contour in 
the lower part of the figure) that corresponds to vanishing determinant of the matrix ( 2 . 11 ): 
— 4:XhXs)v‘^ = 0. In fact it is interesting to consider the special case of — AXhXs = 0, 
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Figure 1. Contour plots for masses of the non-standard (/ 12 ) Higgs particle in the plane {Xh, k) for 
fixed value As = 0.1. The observed Higgs mass is set to be = 125.7 GeV, the vev of the Higgs 
doublet is fixed at u = 246.22 GeV. In the bottom part of the plot {Xh < Xsm = = 0.13) 

the heavier Higgs is the currently observed one, while in the upper part (Xh > Xsm) the lighter 
state is the observed one. The horizontal black line Xh = Xsm separates the two scenarios. White 
regions in the upper and lower parts are disallowed by the positivity conditions for and , 
respectively. 


then if one requires that the potential has a minima for u, 7 ^ 0 it is necessary to assume 
that which implies that the potential could be written as 

- (a^/VP + Af |5p) + (a^/VP + Af |5p)' (2.19) 

with /X]^'^. This potential has equipotential contours (in the unitary gauge) 

on ellipses such that X]^‘^\H\‘^ + A^^ISp = constant. Fluctuation along the ellipses that 
corresponds to the minimum is the massless mode. Note that this is a mode that exists 
even if Ux 7 ^ 0 , so this parameter choice is different from the case discussed below at the 
end of this section where Vx is approaching 0 (so det(Al^) —)■ 0 as well). 

2 

The behavior of Vx is presented in fig 2. For Xh < Xsm and Xh > (as required by 
the scalar mass positivity) one finds 

0<vl<v^^^ ( 2 . 20 ) 

As 
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Similarly for Xjj > XsM and Xh > ^^/(dAg) + Xsm Vx is limited by 


! AgM 

As 


< Vx < OO 


( 2 . 21 ) 


and it is diverging at the parabola Xh = / (4A5) + Xsm- For Xh < Xsm only the physical 

region corresponding to > 0 ( 4 Ai^AsM — > 0 ) is shown in fig 2 . 



Figure 2. Contour plots for the vacuum expectation value of the extra scalar (left 

panel) and of the mixing angle a (right panel) in the plane (Ajy, n) for fixed value Ag = 0 . 1 . In the 
bottom part of the plots the positivity of limits the allowed area by Xh > k^/(4As), so the 
white region is disallowed. The observed Higgs mass is set to be = 125.7 GeV, the vev of the 
Higgs doublet is fixed at u = 246.22 GeV. In the bottom part of the plots {Xh < Xsm) the heavier 
Higgs is the currently observed one, while in the upper part {Xh > Xsm) the lighter state is the 
observed one. The horizontal black line Xh = Xsm = 0.13 separates the two scenarios. 

In order to understand the behaviour of the mixing angle a one can directly adopt the 
formulae (2.13). Since a varies in the range [—7r/2,7r/2], the absolute value of |q;| and its 
sign can be read from the inverse of cos(2a) and sin(2Q;), respectively. The coupling of hi 
eigenstate (the one that has 125.7 GeV mass) to VV is proportional to cos a therefore the 
LHC data favours regions of a ~ 0. As can be found from (2.13) small a corresponds to 
either to t 00 (for Xh > Xsm) or to —>■ 0 (for Xh < Xsm)- One can see the same 
behaviour from fig. 2. 

It is worth investigating the SM limit of the VDM model. Fig. 2 is a good starting point 

as it is easy to recognize regions in the parameter space that imply vanishing corrections 

(relative to the SM) to hi couplings (this is what we define by the SM limit). So, for 

2 

Ajy > Xsm this is the parabola Xh = + Xsm where Vx ^ 00 and a vanishes, while for 

Ajy < Xsm it is the vicinity of Xh = Xsm- The case Xh > Xsm is less interesting as it is 
just the decoupling limit with ~ 2 A 5 r;^ + •••—;■ 00 , see fig. 1, when the effective low 
energy theory is just the SM with h 2 being integrated out. The mixing angle a behaves in 
that region as a ^ k/X s{v/Vx)‘^ + ••• ■ More interesting is the region with Xh approaching 
Xsm from below, as there the mass of /i 2 goes to zero (also —)• 0), as seen in fig. 2. In 
this region the model contains a SM-like scalar with the mass of 125.7 and almost massless 
state (/ 12 ) that decouples from VV, hihi and fermions, however its cubic {Khih^) and 
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quartic scalar couplings remain. This limit of the VDM model turns out to be 

phenomenologically unattractive since the Higgs boson hi would decay invisibly into pairs 


3 The renormalisation gronp eqnations 


As will be illustrated shortly, in order to investigate vacuum stability in the VDM 
model it is necessary to use 2-loop RGBs. We have adopted SARAH [21-24] to obtain the 
full set of 2-loop beta functions. However, in order to facilitate further discussion we show 
below the beta functions of the gauge couplings at 1-loop level, even though whenever we 
impose any constraints on the parameter space that rely on the RGE running we always 
use 2-loop beta functions. 

(3.1) 


Pga = 


we have 


/’ll = 




dV 


PS = -y52 


and = -7gl 


The various Higgs quartic couplings at 1-loop run as 


(3.2) 





27 A 9oo 9^ 90 ^ 9 \ ^ t \ 9 9 

200 ^^ + '^3i92 + g 52 “ 5^1 ~ ^92^h + 24 Aj ^ -|- k 

- 6Tr(v.yiy„yi) - 2Tt(y,y,^y,y^^ - 6Tr(yrfy]yrfy]) 

+ 12AHTr(ydyj) +4AHTr(yeyi) + 12Aj^Tr(y„yj) 
2 ( “ ^^9x^S + 27^^ -|- 40A| -|- 
^ - 9 ( 7 ? - 905 ? - 455 ? + 120 Ah + 8 OA 5 + AOk 

+ 60Tr(yrfy]) +20Tv(yey;) +6014 (y,yj) 


(3.3) 

(3.4) 


(3.5) 


Above denote the corresponding Yukawa matrices. Since here we are mainly con¬ 

cerned with the case of the masses for the extra scalar /12 and Z' of the order of the 
electroweak scale and since initial conditions for RGE running will be specified at Q = mt 
(relatively large scale) therefore we will adopt the above beta functions neglecting decou¬ 
pling of extra degrees of freedom below the scale Q = , Mz' ■ 

We have verified that, in the SM limit, our 2-loop running of \h{Q) agrees with known 
results [34]. 

In order to explore the stability and positivity of the theory we used the two-loop 
RGEs and the tree-level potential. To improve the precision of this work further one would 
likely need the one-loop improved effective potential before extending to 3-loop RGEs. 

An example of a representative point in the parameter space at Q = mt is pictured 
in figure 3, where we show both 1- and 2-loop running of the gauge and scalar quartic 
couplings. A few comments are here in order. The running of gauge couplings (left panel 
of 3) is rather stable and similar for 1- and 2-loop beta functions. However the Yukawa 
coupling Yt already shows (left panel) some sensitivity to the approximation adopted for 
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t=Log,„(Q[GeV]) t = Log,„(Q [GeV]) 


Figure 3. Running of various parameters at 1- and 2-loop, in solid and dashed lines respectively. 
For this choice of parameters Xh{Q) > 0 at 2-loop (right panel blue) but not at 1-loop. Xs{Q) 
is always positive (right panel red), running of k{Q) is very limited, however the third positivity 
condition K{Q) + 2^yXHiQ)Xs{Q) > 0 is violated at higher scales even at 2-loops (right panel green). 


the beta functions. It is important to note the relevance of 2-loop running of Xh{Q)] 
as seen in the right panel of fig. 3 the scale of instability (i.e. a scale Q* from which on 
Xh{Q) is negative) is very sensitive to the RGE running precision. For 1-loop beta functions 
Q* ~ 10^^ GeV while for 2-loop approximation Xh{Q) remain positive up to scale 10^® GeV. 
Other quartic couplings {k{Q) and Xs{Q)) do not require 2-loop beta functions as their 
evolution is nearly the same for 1- and 2-loop beta functions. It is worth noticing that I3\g 
is always positive, so the stability condition > 0 can never be violated radiatively. The 
evolution of k is rather mild as (5^ (x k therefore at least for small k the evolution is quite 
suppressed. In the right panel of fig. 3 we also show the running of the positivity condition 
k{Q) + ‘1\/^h{Q)^s{Q) > 0. Here the 2-loop effects are again important. The 1-loop 
curve terminates already around Q ~ 10^^ GeV, the scale at which 1-loop Xh{Q) becomes 
negative, so that the positivity condition k,{Q) -|- ‘2^\/Xh{Q)Xs{Q) > 0 can not be verified. 
On the other hand the 2-loop running of this condition shows that it is satisfied up to 
Q ~ 10^^ GeV. The choice of initial conditions for this plot illustrates the fact that there 
exist initial conditions (i.e. points in the parameter space (Ajy(mt), K{mt), Xs{mt), gx{mt))) 
such that even though the condition k{Q) + 2y^Ajy(Q)A5(Q) > 0 is satisfied at low scale 
it fails at high energies. 

3.1 Stability 

The constraints (2.7) can be used to determine areas of parameter space 
{XH{rnt), K{mt), Xs{mt), gx{mt)) in which the conditions for stability/positivity of the po¬ 
tential are satisfied at all renormalization scales. In the SM, the absolute stability is ensured 
just by the positivity of the quartic coupling at all energy scales Q: Xh{Q) > 0. Since 
the mass of the Higgs boson is known experimentally the initial condition for running of 
Xh{Q) is fixed as Ajy(mt) = M|^/(2u^) = Xsm = 0.13 and for this initial value Xh{Q) 
becomes negative at some scale causing the instability. However here, in the presence of 
the extra scalar S this is not necessarily the case; the LHG Higgs mass measurement fixes 
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the following combination of couplings and vev’s: 


+ Xsvl ± ^ A|u^ - 2\hXsv‘^vI + \]jv^ + (3.6) 

It is easy to see that the VDM model has the freedom to increase the value of Xh at 
low scales; a freedom which the SM does not possess. Larger initial values of Xh such 
that XH{mt) > XsM are allowed delaying the instability (by shifting up the scale at which 
Xh{Q) < 0 ). There is also another remedy for the instability within the VDM model; even 
if the initial Xh is smaller than its SM value, A//(mt) < Xsm-, still there is a chance to 
lift the instability scale if appropriate initial value of the portal coupling K{mt) is chosen. 
This effect is caused by the positive k? contribution to the beta function that partially 
compensates the negative top-quark effect, see (3.3). Fig. 4 illustrates the way the SM 
stability problem encoded by Xh{Q) < 0 for Q > Q* could be relaxed within the VDM 
model. The white region above the horizontal line XH{mt) = Xsm shows the region of 
XniiT^t) > Xsm so the positivity of Xh up to the Planck scale could be easily guaranteed. 
On the other hand the white region below the line XH{mt) = Xsm shows those pairs 
of {Xh,k) for which even though the starting point for Xh evolution is lower than for 
the SM, nevertheless the extra positive contribution to f3x^ makes Xh positive up to the 
Planck scale. Clearly for large k the stability region increases (for negative k the other 
stability condition gives tighter constraint). The colorful regions show the scale at which 
Xh becomes negative. The three panels shown in fig. 4 correspond to three different pairs 
of initial values for {gx{mt), Xs{mt)). As seen, the sensitivity to those choices is very weak 
even though gx{mt) and Xs{mt) vary in a wide range, in fact this is understandable since 
the evolution of Xh{Q) is influenced by gx{rnt) and Xs{mt) only indirectly through the 
presence of k? in the beta function /3 a^, see (3.3). 


5x(^f)=0 05. As(mf)=0.05 
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Figure 4. The stability frontier for the H direction: these plots identify the renormalisation 
scale t* = Logic (Q*) at which Xh{Q*) = 0 and the vacuum becomes unstable, as a function 
of {X{mt), K{mt)) ■ The horizontal solid black line corresponds to XH{mt) = Xsm — 0.13. For 
t = LogioiQ*) > Logio{Mpi) = 19.09, the vacuum is absolutely stable up to that scale. At lower 
values the vacuum is meta-stable however a lower scale does not imply increased instability, one 
must further determine the tunneling rate. 

In the VDM model the SM stability problem (positivity of Xh) is easily solved as was 
illustrated above. However in this case positivity requires two extra constraints: A 5 > 0 
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and K + 2\/XhXs > 0. Since Pxg > 0 therefore whenever Xs{mt) > 0 the positivity 
is preserved during the evolution. However the second extra condition is non-trivial, as 
illustrated in fig. 3 it is possible that k{Q) + ‘2'x/X h{Q)Xs{Q) changes sign while running 
from low energies up. Fig. 5 shows the scale at which k{Q) + 2^y XHiQ)Xs{Q) becomes 
negative as a function of {X{mt), K{mt)) for three fixed sets of {gx{'nit)-,^s{'fnt))- 



Figure 5. The “in between” stability frontier : these plots identify the scale = Logi[){Q*) at 
which the positivity condition k{Q) + 2yJX h{Q)Xs{Q) fails and the vacuum becomes unstable, as 
a function of {X{mt), K{mt)) for fixed choices of {gx(jnt),Xs(jnt)) specified above each panel. The 
horizontal solid black line corresponds to XH{rnt) = Xsm — 0.13. The gray area is excluded by the 
requirement that there is no Landau poles up to the Planck mass. 

Stability of the U{l)x dark matter model was also discussed in [35]. The vacuum 
stability induced by the dark matter has also been considered in the context of complex [36] 
and real [37, 38] extra scalars serving as dark matter candidates. 

3.2 The Landau poles 

As we have discussed above, the additional freedom in the Higgs sector, that is due 
to the presence of the Higgs portal allows one to increase the low scale value of 

Xh sufficiently to avoid its negative value (instability in the H direction) at high scales. 
However, this possibility is bounded from above by the requirement that there are no 
Landau poles in the evolution of Xh (or any other parameter - a pole in the evolution of 
any coupling implies divergence of all of them at the same energy) up to a chosen high 
scale, e.g. the Planck mass. In fig. 6 we show contour plots of Xh{Mpi) in {XH{rnt), 
space for fixed gx{mt) and Xs{mt). It is clear that too large K{mt) or Aiy(mt) implies early 
divergence of Xh{Q)- Also when gxijnt) and/or Xs{mt) grow (from left to the right panel) 
the safe region shrinks in agreement with expectations. 

A few comments are in order here. Quartic scalar couplings Xh, Xg, k and gx are free 
parameters in the model. However for perturbative expansion to be valid their values can 
not be too large since otherwise the expansion has no chance to converge. The maximum 
adopted by various authors is to some extent subjective and usually varies between 1 and 
d-TT. Here we have the advantage of knowing both 1- and 2-loop beta functions for the RGE 
evolution of the couplings, therefore the relevance of the 2-loop could be quantitatively 
estimated for different (large) values of couplings. We found that requiring the 2-loop 
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Figure 6. Contour plots of Xh{Mpi) in the plane of {X{mt) , K{mt)) for fixed gximt) and As(mt) 
specified above each panel. The horizontal solid black line corresponds to XH{mt) = Xsm — 0.13. 
The plots allow one to identify regions (white) in which the Xh{Q) Landau pole is below the Planck 
scale. 

correction to be smaller than 100% of the 1-loop result for any quartic scalar coupling, i.e. 

< 1, implies that the coupling should not exceed a value close to 27r. 
Therefore, in the numerical results, whenever it is stated that a Landau pole appears, it is 
meant that the corresponding coupling reaches a value of 27r. For larger couplings the 2- 
loop contributions start to dominate, so that one can not trust the perturbative expansion, 
truncated to this order. 




Figure 7. The regions of allowed parameter space (inside the “boomerang” contours) where cou¬ 
plings remain perturbative, positivity is satisfied up to Mpi in the plane (| sina\, Mh^), at 2 loops. 
Left and right panels corresponds to positive and negative k, respectively. The allowed area de¬ 
creases as the magnitude of n increases. The blue area denotes region where hi coupling to vector 
bosons is away from its SM value by more than 15%. The region over dashed red line is excluded 
at 95% CL by the analysis of the Peskin-Takeuchi S,T parameters. 
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The VDM model without the additional U{l)x reduces to a Higgs portal extension of 
the Standard Model, which has been explored recently in [39-41], These papers note that it 
can be useful to consider a change of parameterisation from {Xh, h, Xs) —>• k, sin a), 

in which one can re-express the quartic couplings as follows: 



(3.7) 



(3.8) 


One can then identify the regions in this parameter space, (| sin a], M^^), that are absolutely 
stable (i.e. the stability conditions are satisfied up to the Planck scale) and do not have 
Landau poles at any renormalisation scale before the Planck scale, as pictured in fig. 7. 


We show here those plots in order to compare our results with those of [41], which differ 
in so far as our study include 2-loop (not 1-loop) running, and where our model includes 


the additional parameter gx- 

4 Constrained parameter space 

In this section we show results of scans over the parameter space (Xb, k, Xs,gx}- Usu¬ 
ally Xs will be fixed and specified, gx is either fixed (in figs. 9 - 11) or varied 0.1 < < 1 

while Xff and k will be scanned such that 0 < Xh <0.25 and —0.5 < k < 0.5 (with some 
exceptions when — 1 < «: < 1). If RGE running is employed then the scan range should be 
regarded as for initial values of running quantities at Q = mt- 

4.1 Collider constraints 

In addition to the theoretical requirements discussed above we are going to impose 
some experimental constraints. First of all there exist limits on branching ratio for invisible 
Higgs boson decays. Searches for Higgs decaying invisibly has been carried out by both the 
ATLAS and CMS collaborations at the LHC for various production and decay channels. 
ATLAS [42] considered an associated Higgs production with a vector boson iV = 
or Z), assuming SM production they found an upper limit of 29% at 95% confidence 
level on the branching ratio of Higgs bosons decaying to invisible particles. A search for 
invisible decays of Higgs bosons was also performed by CMS [43, 44]. Assuming Standard 
Model Higgs boson cross sections and acceptances, the observed upper limit on the invisible 
branching fraction was found to be 57% at 95% confidence level. It turns out that within 
the VDM model unless the U{l)x coupling constant gx is tiny or 2Mz' > 125.7 CeV decays 
of the observed Higgs boson hi into Z'Z' would dominate with branching ratio exceeding 
the experimental limits. Similar comments apply for decays of hi into pairs of /12 (for 
the scenario with Xh < Xsm)- Therefore in our analysis we simply exclude points in the 
parameter space such that hi decays to Z'Z' or /i 2 fi 2 are kinematically allowed. 

If Xh < Xsm Higgs boson /12 is light so that it could have been prodeced at LEP. Then 
the LEP limits for e"*“e“ —?• Zh 2 apply and should be imposed. Here we adopt the data 
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from [45] where limits on ZZh\ coupling normalized to the SM ZZh coupling = cos a) 
are tabularised as a function of Mh^. For fixed A 5 the limit on a could be easily translated 
into allowed region in the (k, Xu) plane. 

Higgs boson couplings are being measured at the LHC. The ATLAS [46] and CMS 
[47] collaborations limit e.g. ratios ny of Higgs boson VV couplings normalized to their 
SM values. The conclusion is that the observed Higgs has SM-like couplings to the vector 
bosons. Here we will assume, somehow arbitrarily, that the ratio is limited by 0.85 < 
Ky < 1. Note that because of the orthogonal mixing Ky = cos a can not exceed 1 within 
the VDM model. As we have already mentioned above the SM-like nature of hi favours 
regions of large and small Vx for Xu > Xsm and Xu < Xsm, respectively. The constraint 
that originates from 0.85 < Ky < 1 could be expressed as an allowed region in the ( k , Xh) 
plane both for Xu < Xsm and Xu > Xsm case. 

In order to estimate relevance of electroweak precision data we adopt the Peskin- 
Takeuchi S, T and U parameters [48]^. At the 1-loop level, beyond the SM radiative 
corrections to the SM vector boson progators 5Iiyy are not affected by the presence of the 
Z' boson, therefore 6Iiyy are the same as those found in the analysis of the plain Higgs 
portal [41]. At this order shifts in dliz^ and vanish, therefore the S and T parameters 
can be expressed as 


167r cos^ „ An f 6Uwwi0) (5nz2(0)\ 


(4.1) 


whereas the parameter U is too small to be relevant. Using the fit obtained in [49] we 
found that at 95% confidence level the S and T parameters do not constrain further the 
parameter space. These bounds are entirely embedded in region, where the scalar mixing 
angle is too large or couplings are non-perturbative (see fig. 7). It is worth to emphasize 
that they constrain |sina| only moderately weaker than the full set of the electroweak 
precision observables [41]. In particular, in the important mass range Mh^ > 200 GeV, 
the allowed value of |sina| in S and T approximation is larger by 25%. For low Mh^ 
mass differences grow however that region of parameter space is anyway disfavoured by the 
requirement of absolute stability and excluded by limits on Ky. 


4.2 Dark Matter abundance 

Before we proceed to constrain the parameter space by measurements of DM abun¬ 
dance, in fig. 8 we show results for QuMh'^ as a function of the DM mass {Mz') obtained 
varying coupling constant Qx and choosing a few representative values of other parameters. 

We have calculated Q^Mh^ adopting standard textbook methods for cold relics, see 
e.g. [50]. Relevant vertices, Feynman diagrams and corresponding contributions to Z'Z' 
annihilation cross section for various final states are collected in appendix A. We have 
checked our results for QuMh'^ against calculations done adopting the micrOMEGAsS [31]. 

^Since the scale of new physics in our case is not always much above the electroweak scale therefore the 
S, T and U parameters should be used with some extra care. In particular, for Mh2 < Mh^ they do not 
provide a viable estimation of radiative corrections. Luckily that region of parameter space is not allowed 
by other constraints. 
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It turned out that except resonance regions (such that 2Mz' ~ a mass of s-channel res¬ 
onance) and vicinities of thresholds for new final states, determined via the cold 

relics technique agrees pretty well with the result provided by the micrOMEGAs. How¬ 
ever, in order to have also those special regions under control we have decided to adopt the 
micrOMEGAs hereafter. For illustration, in 8 we compare obtained by adopting 

results contained in appendix A with those from micrOMEGAs. 


Ah=0.127, k =- 0 . 05 . As=0.1 (Mh,=73.1 GeV, v,,=170 GeV) 



A„=0.127, k = 0 . 02 . As=0.1 (^^^=110. GeV, v„=250 GeV) 



A„=0.18, k=0.15, As=0.2 (M„^=205. GeV, v„=300 GeV) 



Mz' [GeV] 

A„=0.145, /(=-0.1, As=0.2 GeV, v„=500 GeV) 



Figure 8. The left and right panels show O.oMh’^ as a function of Mz' for a case of Mh^ < 
125.7 GeV {Xh < Xsm) and > 125.7 GeV {Xh > Xsm), respectively. Masses of Z' that 
correspond to a resonance or an opening of a new final states are shown above each panel. The 
dotted and solid lines correspond to results obtained from the micrOMEGAs and results of cold 
dark matter approximation from the appendix A, respectively. The thick horizontal green lines 
correspond to the observed 5cr result for as in (4.2). 


There is a comment here in order. In fig. 8 and similar that will follow, one sees 
that for large Mz' it is typical that Q.DMh'^ decreases as Mz' grows. In fact it is easy 
to understand such behaviour. In those plots potential parameters are fixed, so is Vx, 
therefore increasing Mz' implies growing gx, so that {av) increases and therefore 
decreases. This fact can be seen in the easiest way by looking at the contribution from 
direct DM-scalar interaction coming from the vertex V^j' oc (appendix A); then (av) oc 
WiWMl oc Qxlv"^ and therefore cc I/Qx- 

4.3 Allowed parameter space 

All constraints that we have considered are collected in figs. 9, 10 and 11. Choosing 
randomly points in the parameter space in the region (white) allowed by perturbativity, 
stability, LEP and LHC data we require that the DM abundance remains within the 5 ct 
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limit of [51] 


VLDuh^ = 0.1199 ±0.0022 


(4.2) 


Points that fit into the allowed region of (4.2) are shown in figs. 9-11 as dots within the 
white region, for them all the constraints are satisfied. 
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Figure 9. Combined plots of allowed and disallowed parameter space in the plane (Ai/(TOt), K{mt)) 
for gxirrit) = gi{rnt) and As(mt) = Xsiyiimt) = 0.13. The thin red line denotes the frontier above 
which a Landau pole of Xh appears below the Planck scale. The thin blue line denotes the absolute 
stability frontier. Below the thin green line the positivity condition fails at some renormalisation 
scale (its wavy shape is a numerical artifact). The yellow region is disallowed by invisible Higgs 
decays. The green area denotes LEP exclusions on Higgs-like scalars. In the outer red area positivity 
fails at the low scale, while in the orange area no physical solution of the vev Vx exists. The blue 
area denotes an excess of the hi Higgs couplings to vector bosons (kv). The remaining allowed 
region is in white. The points (in the white region) coloured with respect to the mass of the extra 
Higgs boson Mh^ are those for which also ^oAth'^ constraint (4.2) is fulfilled. 
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Figure 10. Same as in fig. 9 however colouring of allowed points is here with respect to sina. 


4.4 Direct detection of Dark Matter 

Interactions of VDM with nucleons are mediated by the Higgs particles. For the elastic 
scattering cross section we obtain 

^Z’N = ^gWhNN sin" 2a ( i i J , (4.3) 

\ “hi “h2 / 

where /r = Mz'Mn/{Mn+Mz') is nucleon-DM reduced mass, ghNN = l/u(iV| ^qQQ\^) 
is the effective Higgs to nucleon coupling [52], 

In this subsection we are going to show results for az 'at as a function of Mz' for those 
points (green points in figs. 9 - 11) in the parameter space which satisfy all the constraints 
including the DM abundance. In the following subsection we are showing results obtained 
with micrOMEGAs for A 5 fixed at A 5 = 0.2. 

4.4.1 Light dark matter 

We start with the case of Xh < Xsm- Then, since Vx is limited from above as in 
(2.20) and the scanning interval for Qx is 0.1 < ( 7 ^; < 1, therefore the DM mass is bounded 
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Figure 11. Same as in figs. 10, however for gx{mt) = 0.25 and Xs{mt) = 0.05 (left panel) and for 
gxirrit) = 0.6 and Xsinit) = 0.2 (right panel). 

from above by Mz' < v{Xsm /~ 200 GeV. It turns out that in practice there is no 
consistent points in the parameter space with Mz' > 120 GeV. Note also that in order to 
prevent invisible hi decays hi —)• Z'Z' it is required that Mz' > 62.9 GeV, so the allowed 
range for DM mass is 62.9 GeV < Mz' < 120 GeV. In figs. 12 and 13 we are showing ctz'n 
coloured with respect gx, Mh^ and sin^(2a) in order to learn properties of the points that 
are plotted. We also show experimental limits for (Jz'n from XENONIOO, LUX (2013) and 
anticipated results for XENON IT. As it is seen from the figures, points that are consistent 
with the present data correspond to medium gauge coupling gx ~ 0.5, /12 slightly lighter 
than the observed Higgs ~ 110 d-125 GeV and small mixing angle 0 a ^ vr/S. The 
DM mass varies between 60 and 120 GeV with heavier states favoured. 

Fig. 14 is helpful in order to understand why there is no points in the parameter space 
with Mz' ^ 120 GeV. As it is seen in the right panel of the figure for large Mz' the DM 
abundance Q^uh'^ is very small since almost all annihilation channels (except hihi and 
tt) are open for Z' of that mass (so that the annihilation cross section is large). It is also 
instructive to look at correlations between Mh^ and Mz' shown in the left panel of fig. 14. 
It turns out that there are two regions consistent with the constraints: i) Mz' ~ 65 GeV 
and a) Mz' ~ — 5 GeV. The first one corresponds to the vicinity of hi resonance 

in the right plot in the figure. The dip at 62.9 GeV is quite steep such that deviation 
by about 5 GeV is sufficient to reach O.DMh'^ ~ 0.1. The other side of the resonance is 
excluded by the requirement of no invisible hi decays. In fact it also easy to understand 
the second region that corresponds to the other side of the summit seen to the right of 
the hi resonance in the right panel. There the sudden drop of QoMh'^ is caused by the 
opening of /i 2 /i 2 final state. In a real case, since annihilating Z' pairs are not exactly at 
rest, therefore the /i 2 /i 2 channel opens by 5 GeV earlier. 
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Figure 12. The figure shows the DM-nucleon cross section, az'N, as a function of the DM mass 
Mz' for points which satisfy all other constraints for Xh < Xsm- The singlet quartic coupling is 
fixed at As = 0.2. Colouring corresponds to the strength of the gauge coupling g^- The nearly 
horizontal lines are the experimental limits for az'N from XENONIOO, LUX (2103) and anticipated 
results for XENON IT. 
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Figure 13. Same as in fig. 12, but colouring is with respect to Mh^ and sin^(2a) for the left and 
the right panel, respectively. 


4.4.2 Heavy dark matter 

In this section we are considering the case Xh > )^SM- Then, since Vx is limited from 
below as in (2.21) and the scanning interval for Qx starts at 0.1 therefore the DM mass 
is bounded from below by Mz' > 0.1v{Xsm /~ 20 GeV. However since invisible 
Higgs decays should be prevented therefore Z' must be even heavier, so that here Mz' > 
62.9 GeV. Results for az'N are presented in a similar manner as in the case of the light 
DM, so in figs. 15 and 16 we are showing the cross section plotted against the DM mass Mz' 
with colouring corresponding to Qx, and sin^(2Q;). As it is seen from the figures there 
exist points that lay below the LUX 2013 upper limit, they correspond to medium gauge 
coupling gx ~ 0.2-b0.6, wide range of varying from 130 GeV up to 1000 GeV and small 
mixing angle 0 a tt/8. The allowed DM mass varies between 62.9 GeV to 1000 GeV. 
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Ah = 0.11, K = 0.12, As = 0.2 (Mft, = 75 GeV, i/, = 138 GeV) 



Figure 14. The left panel illustrates correlation between between Mh 2 and Mz', while the right 
one shows predictions for as a function of Mz'- The colouring corresponds to the cross 

section az'N- Above the right box resonances and channels which open as Mz' increases are shown. 
Coordinates in the parameter space (\h,k,\s) and corresponding Mh^ and Vx are shown above 
the right panel. 



Qx 


1.0 

0.8 

0.6 

0.4 

0.2 


0 


Figure 15. The figure shows the DM-nucleon cross section, az'N, as a function of the DM mass 
Mz' for points which satisfy all other constraints for Xh > Xsm- The singlet quartic coupling 
is fixed at As = 0.2. Colouring corresponds to the strength of the gauge coupling gx- The solid 
lines are the experimental limits for az'N from XENONIOO, LUX (2103) and anticipated results 
for XENON IT. 


It is worth understanding the origin of points that compose the blue and green hills 
in the vicinity of 100 GeV in fig. 15. As seen from the figure they correspond to different 
strengths of the gauge coupling constant: the left one is made of points with gx ~ 0.2 
while for the right one the coupling is typically larger gx ~ 0.5. The purpose of right 
panels in figs. 17 and 18 is to illustrate QoMh'^ dependance on Mz' for representative 
points {Xh, k, Xs) in the parameter space. Note in the left panels of the figures that points 
in the space {Mz', Mh^) are grouped into two rough sets, one for Mz' ^ Mh^/2 and the 
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Figure 16. Same as in fig. 15, but colouring is with respect to Mh 2 and sin^(2a) for the left and 
the right panel, respectively. 


other for Mz' ^ M/i2/2. Inspecting the right panel of fig. 17 it is easy to see that the later 
group is composed of points sitting on the left slope of a heavy Higgs (/ 12 ) resonance (at 
M/i 2 = 400 GeV for this particular example). Since the resonance dip is relatively wide 
in this case, therefore the shift above the the nominal resonance is substantial, typically 
of the order of 100 GeV. Then when Mz' increases eventually /i 2 /i 2 channel opens up, 
the cross section increases and drops leading to higher Mz' consistent with the 

observed DM abundance, that explains the other group of points in e.g. fig. 17 located above 
Mz! ~ Mh^/2. The two groups are separated by the presence of the summit in 
between them. The fig. 18 shows that the two groups have different Qx and therefore they 
could be identified as the points that compose the two hills in fig. 15: the blue hill is made 
of points that lay below Mz' ~ Mh^/2 while the green one of those with Mz' ^ Mh^/2. 
The right panel in fig. 18 illustrates the mechanism of the very strong correlation between 
Mz’ and that is observed for large Mz’ along the line Mz' ~ Mh 2 - As it is seen 
from the figure the correlation is caused by the steep drop in that corresponds to 

opening of ^ 2^2 hnal state. 

Constraints coming from the direct detection can be compared with results of [17]. 
However, a detailed comparison is quite complicated, therefore we limit ourself to a con¬ 
clusion that qualitatively results obtained in [17] for the Abelian case agree with those 
found here. The results presented in the fig. 2 therein present similar behaviour to those 
of fig. 17 in this work. It can be seen that for a given Mh^, when Mz' approaches Mh^/^ 
the nucleon scattering cross section diminishes and LUX bounds can be easily satisfied, 
whereas for Mz' between {Mh^ + Mh^)/2 and Mh^ the cross section is substantially larger 
(the points for Mh^ < 600 GeV). Similarly the LUX bounds constrain the vicinity of 
Mz' = 80 GeV. 

5 Summary 

In this paper we explored an extension of the Standard Model gauge symmetry by an 
extra U{l)x factor. The scalar sector of the model consists of a standard SU{2)l doublet 
(H) and a complex scalar (S) charged under this U{l)x in order to provide a mass for 
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Ah = 0.18, K = 0.18, As = 0.2 (M*, = 400 GeV, i/, = 622 GeV) 



Figure 17. The left panel illustrates correlation between between Mh 2 and while the right 

one shows predictions for as a function of Mz>. The colouring corresponds the the coupling 

(Jz'N- Above the right box resonances and channels which open as Mz' increases are shown. 
Coordinates in the parameter space {Xh, k, ^s) and corresponding and Vx are shown above 
the right panel. 



Ah = 0.134, K = -0.054, As = 0.200 (M*, = 700 GeV, v, = 1107 GeV 



Figure 18. The left panel illustrates correlation between between Mh 2 and Mz', while the right 
one shows predictions for as a function of Mz>. The colouring corresponds to the cross 

section gx- Above the right box resonances and channels which open as Mz> increases are shown. 
Coordinates in the parameter space {Xh, k, Xs) and corresponding and Vx are shown above 
the right panel. 


the extra gauge boson. The possible spectrum of scalar states and mixings between them 
have been discussed in detail. Adopting 2-loop RGE running we have required the vacuum 
stability (positivity) conditions to be valid up to the imposed cutoff scale {Q = Mpi) of 
the model. 

We have compared the sensitivity of the running of the scalar quartic couplings and 
the stability of the vacuum on the loop order of the RGEs and found that the 2-loop RGEs 
make a significant effect on the 1-loop result, which could be improved further with higher 
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order corrections. 

We found that an increase in by even a modest amount easily eliminates the 

electroweak vacuum stability problem up to the Planck scale, which is possible due to the 
additional degrees of freedom in this theory. When the Higgs portal is present 

(so /« / 0) then the initial value \H{mt) could be chosen larger than its Standard Model 
value (where it is fixed by the Higgs mass measurement) in this way escaping the dangerous 
possibility of \h{Q) becoming negative. In addition the model discussed here can have a 
stable vacuum even if \H{mt) is chosen below its Standard Model value, since the RGE 
beta function for Xh contains an extra positive contribution proportional to which can 
lift Xh efficiently enough while running up. For this mechanism to work K{mt) must be 
large enough. 

The presence of an extra neutral scalar makes the phenomenology of the scalar sector 
more attractive and richer but still testable at the LHC and future colliders. We have 
focused here on dark matter aspects of the model. The parameter space of the models has 
been investigated in detail and regions consistent with theoretical, collider and cosmological 
constraints have been determined. In particular we have shown that the DM-nucleon cross 
section could be consistent with the LUX and XENONIOO limits, and also looked at why 
this is so. It has also been shown that the anticipated XENON IT limits on az'N may be 
satisfied within this model. 
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A Relic abundance 


In this appendix we present a detailed account of the tree-level calculation of the relic 
abundance in this model. This has been used to compare with the result generated with 
micrOMEGAs. 


A.l Thermally averaged cross-section 

Thermally averaged cross-section is defined as the annihilation cross seciton a times 
Mpller velocity v averaged over the Boltzmann distribution 

f av exp(-A/pf -F M^i/T) exp(-A/p| M^-,/T)(fpi(fp 2 

(av) = - - -. (A.l) 

/ ewi-ypj + MlJT) exp(-Y'pi + M^,/T)d?pid?p2 


In the nonrelativistic limit m ^ T it can be approximated using the formula 

d(s) 


(cj|u|) = 


4M2 




L'(.) - Af) 

2 ^ ^ 4M|, 


- + 0(4 

X 


-I s=iMl, 


(A.2) 


- 22 - 



where x = and cr{s) = 2y^ s{s — 4:M‘^,)a{s) can be written as the sum of contributions 
coming from all possible final states 


d — ^ff ^W+W- + ^ZZ + + ^hih2 + ^h2h2- 

f 


(A.3) 


Cross sections a for fermions and vector bosons (presented in A. 4) can be easily expressed 
as a functions of s only, but for annihilations into scalars the amplitudes squared depend 
also on the Mandelstam variable t. Therefore respective &{s) needs to be written as 
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^hihj (s) — 

he resped 
expressed by s and 9 using 


sin 0 

-/3iMh.,Mh,) de—fis,t)de, 


where f{s,t) is the respective sum of all amplitudes squared for given i,j and t can be 
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+ 2^ 


ls^-4Ml,s^ 

-2s 

(vj + vj/ 

+ 
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+ s2 

1 s \ 

s 


(A.5) 


cos 6. 


In the nonrelativistic limit, using the formula (A.2), this tedious integration of f{s,t) can 
be avoided by changing the order of the limit and the integral [53], then 

1 


dhihAs = 4M^,) = 




(A.6) 


9(1 + 5ij) 

where t{s = 4M|,) = (M|, + M|.)/2 — M^/; note that this expression does not depend on 
COS0, therefore the integration is trivial and gives factor 1. Similarly for the derivative we 
have 

1 

- li( Mu.. M.^.^ 


AihAs — 4M|,) — 


9(1 + <5, 


-(3{Mh„MhA X 


tj) 


+ ((Ml - Mlf + 4(2M|, - Ml - Ml) 


Afis,t) 
dt^ 


J s= 4 M|,, t=t(s= 4 M|,) 


(A.7) 


The presence of the last term comes from the fact, that t^s) = a{s) + b(^s)cos6, where 
a(4M^,) = —1/2 and b{AM‘^,) = oo. On the other hand integral of cos0 vanishes, therefore 
one needs to calculate the second derivative to determine the limit. 


A.2 Relevant vertices 

Scalar mixing matrix 


R = 


cos a — sm a 
sin a cos a 


(A.8) 
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A.3 Useful formulas 

• Lorentz invariant phase space 


1 / 2(mf + m^) {ml-mlY 
/? mi,m2 =^\n- ^ + 2 

OTT V S 




(A.9) 


Sum over polarizations of a vector boson 


N s s 

r,(m) = 3-^ + 


4m^ 


(A.IO) 


Sum over spins of a fermion 
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(A.ll) 


Propagators 
Clip) = 
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(A.12) 
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A.4 Z' annihilation cross section formulae 
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Squares of amplitudes and interference terms that contribute to the 4 diagrams contributing 
to the hihj {i,j = 1,2) final state shown above are listed below where stands for the 
square of the /cth amplitude while fki for {k,l = 1,2,3,4). 

/i(s) = |hif?H(M^0 (A.16) 


/2(s) = |Eh|"H(MzO 
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Ml 


[Ml - Ml - t)l-Ml + Ml + tf 
16M|, 


(A.18) 


\ \Ml+Ml+Ml-s-tf 


2 - 


2{Ml + Ml + 2Ml -s-t) 


Ml 


+ 


{-Ml + s/2f + {Ml + Ml + 2Ml -s-tf + i(M|^ + 3M|, -s-tf + \{-Ml - 3M|, +s + tf 

Ml 


\{Ml +Ml.+ 2Ml -s- t){Ml + 3M|, - s - tf + {s/2 - MI){MI + 3M|, - s - t){-Ml. - 3M|, +s + t) 


2M|, 


iMl + Ml + 2Ml - s - t){-Ml - 3M|, +s + tf {Ml + 3M|, - s - tf{-Ml - 3M|, + s + ^)" 


2M|, 


16M|, 


(A.19) 
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(A.20) 


fi2{s) = 2\Vi^ \EiMz,M^H) 


h3is,t) = 2\Vifvf'vf' 


2t (-M|, + s/2f + ^(M^ - M|, - tf + li-Ml + M|, + tf 
t - M'^, r M|, M|, 


■ Ml, 


(_M|, + s/2)(M2 - M|, - + Ml + t). 


4M|, 


(A.21) 


Ms,t) = 2\Vifvf'vf' 


2 - 


+ Ml + 2M|, - s - i 


M|, 


" ' MI+ Ml+ Ml - s-t\ 

{-Ml + s/2){Ml + 3M|, - s - t){-Ml - 3M|, + s + t) 

(_M|, + s/2)" + + 3M|, -s-t)^ + \{-Ml - ZMl + s + t)\ 

4M|, J 


(A.22) 


/23(s,t) =23?(S;,)|l/fV/' 


- Ml, 


t - Ml Ml 


t (-M|, + s/ 2)2 + i(M/, - Ml - tf + i(-M2 + Ml + t)2 


M|, 


+ 


{-Ml + s/2)(M/^. - Ml - t){-Ml + M|, + t)n 
4M|, J 


(A.23) 


7 ' 7' 1 r ‘2A4 — s — t 

MU.O = 2«(E.)nf tf - -- I —^ ' 


Ml 


{-Ml + s/2){Ml + 3M|, - s - t){-Ml - 3M|, + s + t) 


Ml 


+ 


{-Ml + s/2)2 + \{Ml + 3Ml -s-t)^ + \{-Ml - 3M|, + s + t)"n 


4M|, 


(A.24) 


hA{s,t) = 2\Vl’vff 


2 + 


-Ml - Ml - 2Ml + s 

Ml 


{Ml + Ml + Ml-s-t){t- Ml) 

4(-M|, + s/2)2 + {Ml - Ml - t)^ + {Ml + 3M|, - s - t)^ + {-Ml + Ml + t)^ + (-M?. - 3M|, + s + i)" 


4M|, 


{-Ml + s/2){Ml - Ml - t){-Ml + Ml +t) + {-Ml + s/2){Ml + 3M|, - s - t){-Ml - 3M|, + s + t) 




4(-M|, + s/2)2 + {Ml - Ml - t)^ + {Ml + 3M|, - s - t)^ + \{-Ml + Ml +1)^ + {-Ml - 3M|, + s + 


4M|, 


(A.25) 


B Model files for SARAH 

In this appendix are included the model files used for SARAH [21-24], to study the 
model discussed in this paper. Currently SARAH does not implement a Z 2 symmetry for 
the imaginary scalar component to give S M- S*, nor for the Z 2 for the 17(^1)x gauge field. 
Instead the kinetic mixing couplings have been set to vanish, gi^x = 5a ;,1 = 0, to preserve 
this symmetry. 
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B.l SMVDM.m 


(*-*) 

(* Particle Content *) 

(* - *) 


(* Gauge Superfields *) 


Gauge [[1]] = {B, U[l] 

, hypercharge, 

gl , 

False } 




Gauge[[2]] = {WB, SU[2] 

, left , 



g2, 

True } 




Gauge[[3]] = {G, SU[3] 

, color 

, 


g3, 

False } 




Gauge[[4]] = {X, U[l] 

, darkcharge 

> 

gx , 

False } 




(* Ghlral Superfields 

*) 









FermionFields [ [ 1 ] ] = 

{q, 

3 , 

{uL, 

dL} 

, 

1/6, 

2 , 

3 , 

0}; 

FermionFields [ [ 2 ] ] = 

{1, 

3 , 

{vL , 

eL } 

, 

-1/2, 

2 , 

1 , 

0}; 

FermionFields [[ 3 ] ] = 

{d, 

3 , 

conj 

dR] 

, 

1/3, 

1 , 

-3, 

0}; 

FermionFields [[ 4 ] ] = 

{u, 

3 , 

conj 

uR] 

, 

-2/3, 

1 , 

-3, 

0}; 

FermionFields [[ 5 ] ] = 

{e , 

3 , 

conj 

eR] 


1 , 

1 , 

1 , 

0}; 

(* FermionFields [[ 6 ] ] 

= 

{v, 

3 , conj [ 

vR] 

0, 

1 , 

1 , 

0}; 

ScalarFields [ [ 1] ] = 

{H, 

1 , 

{Hp, 

HO} 

, 1/2, 2, 

1 , 

0}; 


ScalarFields [ [ 2] ] = 

{S, 

1 , 


hs 

, 

0 , 1 , 

1 , 

1}; 



(*- 

(* DEFINITION 

(*- 

NameOfStates={GaugeES , EWSB} ; 

(* - Before EWSB - *) 

DEFINITION [ GaugeES ] [ DiracSpinors] = { 
Fdl -> {dL, 0}, 

Fd2 -> {0, dR}, 

Ful —> {uL, 0}, 

Fu2 -> {0, uR}, 

Fel —> {eL , 0} , 

Fe2 —> {0, eR}, 

FV —> {vL, 0} } ; 


-*) 

*) 

-*) 


DEFINITION [ GaugeES ] [ Additional]= { 

{ LagHC, { Overwrite—>True , AddHG->True}} , 

{LagNoHC, { Overwrite—>True , AddHG->False } } } ; 

LagNoHC = (mu2 conj [H] . H — lambdaH conj [H] . H. conj [H] . H + nu2 conj [ S ] . S 
— lambdas conj[S].S.conj[S].S — L3 conj[S].S.conj[H].H) 

LagHG = -( Yd conj [H] . d . q + Ye conj [H] . e . 1 + Yu H. u . q ) ; 

(* Gauge Sector *) 


DEFINITION [EWSB] [ GaugeSector] = 
{ 


{{VB,VWB[3]} ,{VP,VZ} ,ZZ} , 

{{VX}, {VZp},ZZp}, 

{{VWB[1] ,VWB[2]} ,{VWp, conj [VWp]} ,ZW}}; 











DEFINITION [EWSB ] [ VEVs] = 

{ {HO, (vSM, 1/Sqrt[2]}, (sigmaH , \ [ Imaginaryl ] / Sqrt [2] } , { phiH , 1/Sqrt[2]}}, 

(hs, (vX, 1/Sqrt[2]}, (sigmaS , \ [ Imaginaryl ]/Sqrt [2] },{ phiS , 1/Sqrt[2]}}} 

DEFINITION [EWSB] [ MatterSector] = 

{(* {{vL , conj [vR] } , {VL,ZM}} *) 

{{phiH ,phiS},{hh,ZH}}, 

{{sigmaH , sigmaS } ,{Ah,ZA}} , 

{{{dL}, {conjfdR]}}, {{DL,Vd}, {DR,Ud}}}, 

{{{uL}, {conjfuR]}}, {{UL,Vu}, {UR,Uu}}}, 

{{{eL}, {conjfeR]}}, {{EL,Ve}, {ER,Ue}}}}; 


-*) 

*) 


DEFINITION [EWSB] [ DiracSpinors] = { 


Ed - 

->{ 

DL, 

conj [DR]} , 

Fe - 

->{ 

EL, 

conj [ER]} , 

Fu - 

->{ 

UL, 

conj [UR]} , 

Fv - 

->{ 

vL , 

0}}; 

DEFINITION [EWSB] [GaugeES] = { 

Fdl 

->{ 

FdL, 0}, 

Fd2 

->{ 

0 , 

FdRj , 

Ful 

->{ 

Ful, 0}, 

Fu2 

->{ 

0 , 

Fu2} , 

Eel 

->{ 

Eel, Oj, 

Fe2 

->{ 

0 , 

Fe2}}; 


(*- 

(* Dirac —Spinors 

(*- 


B.2 SPheno.m 


(*-*) 

(* MINPAR *) 

(*-*) 

OnlyLowEnergySPheno = True; 

MINPAR= { {1, lambdaHINPUT } , 

{2 ,L3INPUT} , 

{3,lambdaSINPUT} , 

{10, gxINPUT}, 

{20, vXInput} }; 

ParametersToSolveTadpoles = {nu2,mu2}; 

RealParameters = {vSM,vX}; 

Boundary LowScaleInput={ 

{gx, gxINPUT}, 

{glx,0}, 

{gxl,0} , 

{lambdaH , lambdaHINPUT } , 

{lambdas , lambdaSINPUT} , 

{L3, L3INPUT}, 

{vX, vXInput}}; 


Boundary LowScaleInput={ 






{vSM, Sqrt [4 mz2/(gl *2+g2 “ 2) ] } } ; 

ListDecayParticles = {Fu , Fe , Fd , hh , VZp, Hp} ; 

List Decay P art iclesSB = {{Fu,”Fu. £90 ” } , { Fe , ” Fe . f90 ” } ,{ Fd, ” Fd . f90 ” } } ; 
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